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( ELECTRONIC CIRCUIT SIMULATION ) 



INPUT CIRCUIT DESCRIPTION 
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REGISTER NON-ZERO ELEMENTS CORRESPONDING TO ELECTRONIC 
ELEMENTS OF CIRCUIT DESCRIPTION WITH MATRIX A 



MATRIX REORDERING 
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PRODUCE OPTIMUM T I ME^HrHHMBl At TO DETERMINE 
NEXT TIME POINT t n+1, WHERE t n+1 = t n + At 



S205 



PRODUCE NEWTON EQUATIONS (LINEAR SIMULTANEOUS EQUATIONS) 



APPLY GAUSSIAN ELIMINATION TO SOLVE 
LINEAR SIMULTANEOUS EQUATIONS 



S206 




S2.07 

SOLUTIONS CONVERGED ° m - 
~~~~ YES 



U.S.Appl. Serial No. 09/822,512 Filing Date: " *y ^/c/^^ Applicant Koutaro HACHIYA 

JBk METHOD AND APPARATUS FOR MATRIX REORDERING J^ptECTRONIC CIRCUIT SIMULATION 



FIG. 2 

( START ) 



i:=1 

N:=SIZE OF MATRIX 
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initialized k o^k^N) TrrrrrrTFr^ 



SELECT MINIMAL DEGREE mindeg FROM AMONG 
DEGREES OF PIVOTS RANGING FROM i TO N 



SELECT FROM AMONG PIVOT i TO PIVOT N, 
PIVOT p WHOSE DEGREE IS UNDER mindeg+o? 
AND WHOSE CRITICAL PATH LENGTH IS MINIMAL 



REPLACE PIVOT p WITH PIVOT i 
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WITH RESPECT TO ALL OF ELEMENT e WHOSE VALUE 

yM**£9j&QEBfflm. BY ELIMINATION PROCESS OF PIVOT i, 

CRITICAL PATH LENGTH len(j) OF PIVOT j TO WHICH e 

BELONGS IS fSERgm AS len(j) :=max{len(i)+w(i), len(j)} 
upbATeb 



ADD NON-ZERO ELEMENTS, WHICH ARE PRODUCED BY 
ELIMINATION PROCESS OF PIVOT i, TO COEFFICIENT MATRIX 



i:=i+1 
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for k=1 :n-1 



ELIM INATE 
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MULTIPLICATION 



for i =k+1 : n 

if a(i,k)<>0 then 
a(i.k)=a(i,k)/a(k. k) 
for j=k+1 :n 

if a(k. j)<>0 then 

a(i, j)=a(i, j)-a(i, k)a(k, j) 

end 

end 

end 

end 
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NON-ZERO ELEMENTS 
NON-ZERO ELEMENTS 
(NUMBERS IN CIRCLES REPRESENT 
NUMBERS OF PIVOTS BEFORE 
REORDERING) 
ZERO ELEMENTS 
NON-ZERO ELEMENTS ADDED 
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